Dense loops, supersymmetry, and Goldstone phases in two dimensions 
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Loop models in two dimensions can be related to O(JV) models. The low-temperature dense-loops 
phase of such a model, or of its reformulation using a supergroup as symmetry, can have a Goldstone 
broken-symmetry phase for N < 2. We argue that this phase is generic for —2 < N < 2 when 
crossings of loops are allowed, and distinct from the model of non-crossing dense loops first studied 
by Nienhuis [Phys. Rev. Lett. 49, 1062 (1982)]. Our arguments are supported by our numerical 
results, and by a lattice model solved exactly by Martins et al. [Phys. Rev. Lett. 81, 504 (1998)]. 



It has long been understood that theories of an TV- 
component scalar field <p with O(N) symmetry in Eu- 
clidean space of dimension d can be related to statistical- 
mechanics models in which configurations of loops are 
given Boltzmann weights depending on their lengths and 
intersection properties Q. In the N — > limit, unwanted 
closed loops vanish Q] , and the usual self-avoiding walks 
(polymers) are related to the generic critical point in the 
0(N — ► 0) theory. The N — > limit can be avoided 
by using a model with a Lie superalgebra as symmetry 
( "supersymmetry" ) Q] . Other values of N are also of in- 
terest. In d = 2, exact results based on lattice models ||] 
show that there also exists a massless low-temperature 
(low-T) regime of the model, for —2 < N < 2. The prop- 
erties of this phase — in which, geometrically, the loops 
become dense — are well-understood using Coulomb gas 
techniques pjj. 

In this Letter, we revisit the low-T phase of the 0(N < 
2) models in d — 2. The supersymmetry approach, which 
can be generalized to all integer N, makes it clear that 
for N < 2 there should be a massless phase described by 
the Goldstone modes of the broken O(N) or supersym- 
metry, with very simple scaling dimensions. We argue 
that the usual dense- loops phase || , in which loops never 
intersect, is not this Goldstone phase (in contrast to ear- 
lier expectations j^]): in fact, it is not generic, and an 
arbitrarily-weak perturbation that allows loop intersec- 
tions causes a cross-over to the generic Goldstone phase 
for —2 <N< 2. We present numerical results that sup- 
port our interpretation. There is also a soluble model 
that appears to be in the Goldstone phase || . 

The generic continuum 0(-/V)-invariant action S for the 
scalar field <f> includes the interaction term —X((j)-(j)) 2 (by 
convention the Boltzmann weight is e s ). In terms of 
loops, this is the model introduced by de Gennes EJ, in 
which the loops can cross, but for A > these cross- 
ings are disfavored. It has a second-order transition for 
d > 2, but for d = 2, a transition occurs only for N < 2; 
for N > 2, there is only the high-T or massive unbroken- 
symmetry phase. (For N not a positive integer, the the- 
ory is supposed to be defined by analytic continuation 
of its perturbation expansion. For N integer, this can 



be made rigorous by the use of supersymmetry, as we 
describe below.) When —2 < N < 2, the phase transi- 
tion is second order, and the associated critical exponents 
have been conjectured by Nienhuis H) an d others ||. For 
N < —2, the transition is expected to be first order. 

The existence of the transition for N < 2 implies the 
existence of a low-T phase. From the point of view of 
the generic O(N) <f> A theory, at low T the symmetry is 
broken to 0(N — 1) according to Landau theory, and the 
low-T phase would be expected to be a Goldstone phase 
with massless excitations, described by a nonlinear sigma 
model with target manifold 0(N)/0(N - 1) S^ 1 , a 
sphere. For N — 2, the 0(2) symmetry is not really bro- 
ken, but there is a power-law phase with continuously- 
varying exponents, the Berezinskii-Kosterlitz-Thouless 
(BKT) phase §. For N < 2, use of the perturbative beta 
function in the S N ~ X sigma model, defined by analytic 
continuation in N, implies that the coupling between the 
Goldstone modes in the sigma model becomes weak at 
large length scales ||, like Goldstone phases in higher 
d. This yields a low-T phase with simple forms for the 
scaling dimensions. This is possible in d = 2 for N < 2 
because the Mermin- Wagner theorem no longer applies 
when the continuous SO(iV) symmetry cannot be real- 
ized as unitary operations on a vector field. We discuss 
this theory more fully below, using supersymmetry. 

In contrast, Refs. |^|| defined a particular lattice mo- 
del of strictly non-crossing loops. This was done by trun- 
cating the high-T expansion of a lattice version of the (j) 4 
theory as above ||, in order to simplify calculations of 
critical properties while hopefully remaining in the same 
universality class. The model has Boltzmann weight 

e s = Ij(l + ^-^), (1) 

m 

where i, j are vertices on the honeycomb graph, (ij) de- 
notes an edge of the graph, and the variables Si are N- 
component fixed-length spins, Sf — N . The partition 
function can be evaluated as Z — N L K E where the 
sum is over graphs consisting of self-avoiding mutually- 
avoiding loops, L is the number of loops, E the number 
of edges they occupy, and K ~ 1/T. The critical K for 
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this model is known exactly, K c = (2 + y/2 — N) for 
N < 2 ||. The transition is second-order for -2 < N < 2 
only, and the critical exponents at this critical point, 
which is known as dilute loops (or as self-avoiding walks 
for N — > 0), are known by a variety of techniques |S@ . 
The low-T region for each —2 <N< 2 is attracted to 
a unique massless phase (that varies with N) for all 
K c < K < oo. This phase has nontrivial scaling dimen- 
sions, which are again known PUlcfl. It has become known 
as dense loops (or as dense polymers for N — > 0). In the 
limit N — ► 2, this phase coincides with the critical end- 
point of the BKT phase (which appears in the N = 2 
model for K > l£" c ). For A < —2, we expect that the 
dense-loops phase is massive. 

Both the critical and dense- loops phases are believed to 
be universal as continuum phases, in the sense that they 
are independent of the choice of lattice used in defining 
them. Recent work Oj on the T = limits of certain 
lattice models identifies other massless phases of non- 
intersecting fully-packed loops, that depend on the lattice 
used. We are not concerned with these here, nor with the 
transition for N > 2 in the region \K\ > 1/N, where the 
Boltzmann weight (|l|) can be negative |l2] ]. 

From the critical exponents for dilute loops, it is natu- 
ral to hope fl , that they are in the same universality class 
as the transition in the generic O(N) 4> 4 theory. The idea 
is that any weak repulsion of the lines leads to a crossover 
to the infinitely-repulsive limit in which the loops never 
cross, so that the latter fixed point governs the transition 
for all A > 0. We know of no reason to doubt this. As a 
check, we can reintroduce intersections of the lines in the 
dilute critical theory The graphs of the expansion of the 
generic O(N) model that are not present in the truncated 
model would involve in particular crossings of the loops, 
or multiple occupancies of edges. In the continuum de- 
scription, the most relevant (and most generic) of these is 
where two lines cross. This may be included as a pertur- 
bation of the dilute critical theory by adding to the action 
the integral of a multiple of the corresponding scalar op- 
erator, the so-called 4-leg operator, which has conformal 
weight H hi = g — (g — l) 2 /(4g), where we parameter- 
ized N — — 2 cos7rg, g e [1, 2). We see that this operator 
is marginal at N = 2 and irrelevant (i.e. hi > 1) for 
—2 < N < 2. This strongly suggests that allowing more 
general interactions or configurations of loops in the lat- 
tice model would not change the universality class, at 
least when the perturbations are small, and hence that 
the dilute critical point is indeed generic. 

The model of non-intersecting loops was obtained by 
truncating the high-T expansion, and this would be ex- 
pected not to change the behavior for sufficiently high T 
(small K); apparently this remains true down to T = T c . 
However, this leaves open the possibility that the low-T 
behavior of the models is different, which we will now 
examine. The dense-loops phase is described by a theory 



similar to that of the dilute case, with the same parame- 
terization of N and h± as above, but with g E [0, 1]. One 
sees that the 4-leg operator is now relevant in the whole 
region -2 < N < 2. As N \ -2, hi \ -oo. Since the 
4-leg perturbation is relevant in the low-T phase, it is 
dangerously irrelevant in the critical theory. 

It is thus clear that the dense-loops phase is non- 
generic from an 0(iV)-inodel point of view. In d > 2, 
de Gennes' model appears generic, as the topological ef- 
fect of strict non-intersection, present in d = 2, is absent 
(though in d = 3, there remains the topology of knots and 
linking instead). Even a model in d > 2 of strict non- 
intersection, when reduced to d = 2 by confining d — 2 
dimensions to finite intervals, gives an effective model in 
which intersections still occur. There is also a symme- 
try interpretation for the non-genericity: the dense phase 
of strictly non-intersecting loops has a larger symmetry, 
U(A) instead of O(N), which is based on the possibil- 
ity of consistently orienting all the loops in the partition 
sum [^[[l . Though this symmetry may be broken by the 
boundary conditions, the massless dense-loops phase ap- 
pears to be most easily understood in terms of a non- 
linear sigma model with this larger symmetry ||. The 
4-leg operator is a symmetry-breaking perturbation, that 
reduces the symmetry to the generic O(N). It would be 
surprising if the large-length-scale effect of this relevant 
symmetry-breaking perturbation led back to the same 
higher-symmetry phase. We note that the importance of 
the self-crossings has been observed in the related context 
of Lorentz lattice gases [[l3| . 

We now argue that the flow induced by the 4-leg per- 
turbation leads to the generic Goldstone phase of the 
(j) 4 theory, and begin by examining the structure of this 
phase, before giving numerical evidence that such a flow 
does in fact occur. 

The loop (high-T) expansion of the O(N) model can 
be reproduced by considering instead a model with 
OSp(m|2n) symmetry, with A = to — 2n\ in particular, 
each closed loop incurs a factor of A Wm. The model 
then makes complete sense even for N non-positive, 
(but N must be an integer). The phase diagram, beta 
functions, and scaling dimensions are the same for all 
to at fixed N (though multiplicities of operators may 
vanish for m small ||). The Goldstone phase should 
thus be described by a sigma model with target space 
OSp(m|2n)/OSp(m - l|2n) = S"™- 1 ! 2 ™, a supersphere, 
and a single coupling constant g a The perturba- 

tive beta function is, to leading order, = dg a / dhiL oc 
(N — 2)g%., so for N < 2, the model flows to weak coupling 
for g a > (which is the expected physical sign). Again, 
the Mermin- Wagner theorem does not apply here, and 
symmetry-breaking is allowed for N < 2. The weak- 
coupling fixed point is a theory of free massless scalars, 
with to — 1 bosonic and 2n fermionic components, so is 
conformal with central charge c = to— 1— 2n = N—l, as it 
would be for S' Ar_1 . There will be logarithmic corrections 
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due to the marginally- irrelevant coupling g a . We empha- 
size that in the sigma model, OSp(m|2n) symmetry does 
not allow any other marginal or relevant couplings to be 
added to the action, so that this fixed point is robust, 
unlike the dense-loops theory above. 

Explicitly, the target manifold can be parameterized 
by commuting coordinates Xi, i = l,...,m, and anti- 
commuting coordinates rjj, j = 1, . . . , 2n, subject to the 
constraint J27Li x 1 + Y^j=i V2j-iV2j = 1- The action of 
the sigma model is 



S= - 



1 



d 2 r 



i=i j=i 



(2) 



The simplest example is m = n = 1 (N = —1), which 
corresponds to the S ^ 2 supersphere. The constraint can 



be solved as xi 



\r\\r\2- After elimination of x\ and 



some rescaling, one finds 



S 



d 2 , 



(3) 



At long wavelengths the coupling g a renormalizes to- 
wards zero, and the theory becomes free symplectic 
fermions, with central charge c = —2, as claimed. In gen- 
eral, for N < 1, the partition function of the fixed-point 
theory vanishes, Z — 0. Similar behavior is well-known 
in the dense-polymer (N = 0) phase. 

The supersymmetry analysis is valid for N integer. For 
general N, we are forced again to interpolate formally. 
The availability of the supersymmetry results greatly 
strengthens our confidence in this procedure and in the 
existence of the Goldstone phase for all N < 2. 

Both at the dilute critical point and in the dense-loops 
phase, the central charge is c = 1 — Q(g — l) 2 /g, in the 
parameterization above. For dense loops, g € [0, 1], this 
coincides with c = N — 1 for N = 1, 2. The N — 2 case is 
the BKT phase with c — 1 throughout, and not of interest 
here. We claim that, for all —2 < N < 2, the perturbed 
dense-loops phase flows to the Goldstone phase, which 
is a distinct massless phase even when it has the same 
c. For 1 < N < 2, c decreases during the flow, but for 
—2 < N < 1 it increases. This latter behavior would 
not violate Zamolodchikov's c-theorem, as the theories 
involved are non-unitary. In particular, the N = (c = 
—2) dense-polymer theory should flow to a c = —1 free 
theory. Note that for N = 1 the partition function is 
Z = 1 (after removing nonuniversal constants) in both 
the dense- loops and Goldstone phases, a consequence of 
c = 0. 

Some particularly interesting operators in loop models 
in general are the fc-leg operators (also known as fuseaux 
or watermelon operators), insertions of which give the 
probability that k lines terminate at the same point. In 
the general O(N) point of view, the leading such operator 
would be represented by k insertions of the field <p at the 



same point (with no derivatives), and must necessarily 
be in a totally-symmetric rank-fc tensor representation of 
O(N), which can be assumed to be traceless on all pairs 
of indices. For the OSp(m|2n) formulation, this becomes 
a supersymmetric tensor, that vanishes when contracted 
using the same invariant bilinear form as in the constraint 
in the sigma model. It is useful to choose the index values 
of the set of tensors in a correlation function to occur in 
distinct pairs, to force the lines in the loop model to 
connect the positions of the operators in some specified 
way. This is possible if m is sufficiently large, a typical 
situation in applications of supersymmetry. 

At the dilute critical point and in the dense-loops 
phase, there are nontrivial scaling dimensions Xk (or con- 
formal weights hf. = X]J2) for the fc-leg operators, which 
are known exactly |^JlQ|| (for example, /14 was given ear- 
lier). (In these phases, the symmetric-tensor fc-leg op- 
erators are in general degenerate with others [||], and in 
the dense-loops phase, we assume there is a singlet 4-leg 
operator, which we used as the perturbation above; it is 
the most relevant invariant scalar operator.) In the Gold- 
stone phases, on the other hand, the symmetric tensors 
correspond to functions like spherical harmonics on the 
target space, and their scaling dimensions tend to zero 
as g a — > 0. Since the symmetry is broken for N < 2, the 
precise form of each operator depends on which compo- 
nents are chosen. If the expectation of the field is x\ = 1, 
all other Xi and rjj = 0, then any component of the fc- 
leg operator can be rewritten in terms of the (rescaled, 
free) m—1 x^s and the 2n rjj's, and the correlation func- 
tions will contain various logarithmic factors depending 
on which components are chosen. It is clear that the scal- 
ing dimensions of the leading fc-leg operators are all zero, 
and there will be subleading operators with integer con- 
formal weights (the precise form of correlation functions 
may require more attention to the limit g a — > 0). When 
the dense-loops phase is perturbed by the 4-leg operator, 
it should be possible to see the effective scaling dimen- 
sions of the fc-leg operators (say, in two-point functions) 
cross over to zero, for all fc. 

Numerical transfer matrix calculations were done on a 
square-lattice cylinder of circumference L. The edges of 
the lattice are packed with one line per edge; two lines 
can meet at each vertex in three ways, the weight being 
1 for a non-crossing and w for a crossing. The natural 
0(A^)-invariant boundary condition for a periodic system 
is that the fields <fi are periodic, both for Bose and Fermi 
components. This implies that a loop that winds the 
periodic direction is given the same factor N as a topo- 
logically trivial loop. In a conformally-invariant system, 
c can be extracted from the scaling of the free energy per 
vertex, f(L) = f(oo) — ttc/(6L 2 ). We therefore use this 
formula to extract the effective value of c in a crossover, 
using f(L) at pairs of values of L. Results for N = (see 
Fig. ^) show a clear trend away from c = — 2, and are 
consistent with a slow approach to c = — 1 for L — > 00. 
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FIG. 1. Measures of the central charge vs. w for N = 0. 
The conjectured value for L — > oo is c = — 1 for all w > 0. 

For the fc-leg scaling dimensions, we use the lowest 
transfer matrix eigenvalues Ej~ in each sector with k lines 
propagating along the cylinder. From conformal invari- 
ance jffj] , we define effective values of Xk at finite L from 
Ek — Eq = 2irXk/L; see Fig. || (k = 2 is not included as, 
for N = in our model, Ei = Eq for all L and w). The 
results are consistent with Xk = in the fixed-point the- 
ory. In a situation with a marginally-irrelevant operator, 
one expects Jl5| corrections to scaling to be 0(1/ Lin L), 
which presumably explains the slow convergence. 
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FIG. 2. Measures of the fc-leg exponents Xk = 2hk vs. w 
for N = 0. 

A useful check of our analysis is obtained by modi- 
fying the boundary conditions, in a similar way as in 
Rcf. If we make the Fermi components of the (f> 

fields obey antiperiodic boundary conditions, this will 
give a non-contractible (winding) loop the usual weight 
w c = m — 2n = N if its winding number is even, and 
w Q = rn + 2n = N + An if it is odd. In the Goldstonc 
phase, this will give the n pairs of symplectic fermions 
the antiperiodic boundary condition, while leaving the 
bosonic scalar fields periodic. In the formula for f(L), 
c will be replaced by c — 24/i tw , where h tw is the lowest 
conformal weight in the sector with boundary conditions 
twisted as described. We expect c — 24/i tw = m + n — 1, 
which can be understood as the contribution of the twist 
operator of weight h = —1/8 for each symplectic fcrmion 
pair. Thus we expect 



24/j, 



N - 1 + 3n 



3w Q + W e 



(4) 



for general TV, n. This appears to be confirmed numeri- 
cally for TV — w e = and several values of w Q in Fig. ||. 
Notice how this differs from the usual loop gas behavior 
(w = 0), where c-24/i tw = 1 - 6 arccos 2 (^)/(7r 2 ff ) @. 

For N < —2, we believe that the dense- loops phase is 
massive, and hence flows to the generic Goldstone phase 
only when w is greater than a positive critical value. 
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FIG. 3. Measures of c — 24/i t w vs. w for N 
lines show the prediction of Eq. (0). 
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Finally, it is important that a version of the 
OSp(m|2n)-invariant model on the square lattice is inte- 
grable at w = (2 — N) /4 and can be solved by the Bethe 
ansatz ||. The model is massless for N < 2, and there 
is strong evidence that the central charge is c = JV — 1, 
and that the lowest scaling dimension, other than that 
of the identity operator, is zero; the authors conjecture 
that these are the exact values, but do not identify the 
phase with the Goldstone phase. 

To conclude, we have argued that crossing of loops is 
a relevant perturbation in the dense-loops low-T phase 
for —2 < N < 2, and that the long-distance behavior is 
governed by the Goldstone phase instead, by analytical, 
numerical, and exact approaches. 
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